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Moments of the quark density distribution, moments of the quark helicity distribution, and the tensor charge 
are calculated in full QCD. Calculations of matrix elements of operators from the operator product expansion 
have been performed on 16^ x 32 lattices for Wilson fermions at /3 = 5.6 using configurations from the SESAM 
collaboration and at /3 = 5.5 using configurations from SCRI. One-loop perturbative renormalization corrections 
are included. Selected results are compared with corresponding quenched calculations and with calculations using 
cooled configurations. 



1. INTRODUCTION 

Given the detailed experimental knowledge of 
the light cone distributions of quarks and gluons 
in the nucleon, it is of interest to use lattice QCD 
both to calculate the quark and gluon structure of 
the nucleon from first principles and to reveal the 
underlying mechanisms giving rise to this struc- 
ture. Using the operator product expansion, it is 
possible to calculate moments of quark distribu- 
tions, and we report here the first calculations in 
full QCD §. We also compare full QCD results 
with quenched QCD and with configurations that 
have been cooled to remove all the gluon contri- 
butions except for those of instantons. 

2. DEFINITIONS 

The moments of the spin-independent struc- 
ture functions Fi{x,Q'^), F2{x,Q'^) and the 
spin-dependent structure functions gi{x,Q'^), 
g2{x,Q^) are related to products of Wilson co- 
efficients C„((5^//i^) times hadronic matrix ele- 
ments: 
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J^dxx-F^ix,Q^) = lC:^^{p{x-q){^,) 

dxx" F2{x, 02) = Cl+^{9^) 

S'^dxx-g,{x,Q^) = \C-{^)2{x-Aq){ti) 

J^dxx"g2ix,Q^) = l-2L.[Cf,{^)dni^^) 

-C^{p2{x-Aq)]. 

These matrix elements, denoted {x"q) (which 
equals vlfl_i in the notation of ref.Q), {x"'Aq) 

(which equals ^ai?^ and d, as well as the mo- 
ments of the transversity distribution h(x,Q^), 
{x"'Sq), are related to expectation values of the 
following operators in the proton ground state: 

= i(P5| ^J{^,D^, ■ ■ ■ D^„)i>\PS) 

f(x"A<7)5{.P^,-..P^„} 

= -{PS\ (f)" V;757{.5p, • • • B^„}V^|P5) 
id„S'[^P{^^]---P^,j 

= -{PS\ (1)" V^^g^t^i^^^^] . . . D^^}i^\PS) 
-l-{x-5q)S^,P{,^P,,---P,^} 

= {PS\ (§)" i/^Ts^Ml^^Mi • • • D^„}tl^\PS). 

In the parton model, these matrix elements cor- 
respond to the following moments of quark longi- 
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Figure 1. Contamination by excited states, P, as 
a function of the source RMS radius 



tudinal and transverse distributions: 

= /q dx x" {qi (x) + q^ (x)) 



a;"Ag = Ji^dxx'' {qi{x) - qi;{x)) 
x"-Sq = dxx'' {qT{x) - q±{x)). 

3. SOURCES 

Connected diagrams are calculated using se- 
quential propagators generated by the upper 
two components of the nucleon source J" — 
w>f(C75)/3,/3'df The overlap with the 
physical proton ground state was optimized us- 
ing Wuppertal smearing to minimize the con- 
tamination P = En^o\iJ\^)\^/\{J\0)\^- Figure I 
shows that varying the smearing reduced P by 
over 4 orders of magnitude, yielding an overlap 
with the physical ground state of approximately 
70%. Dirichlet boundary conditions were used for 
quarks in the t-direction. 

4. OPERATORS AND PERTURBATIVE 
RENORMALIZATION 



Table 1 

Operators used to measure moments of quark dis- 
tributions 







mix 


V 


lattice operator 


(a) 

xqc 

u.qc 


6+ 

3 

3+ 


no 
no 


7^ 



QK^ Ddiq 
q jiUAq 

- 3 l-n=lJll>.Diq 


X qc 




yes 




qj{iDiD4}q 
2 2^1=2 97|*^* "*}^ 


x-^qc 




no* 




q^{iDiD4Di}q 

+ qi{2D2D3D^q 
-(3-4) 


Aqc 


< 


no 





97^739 


xAqi"^ 




no 




g7^7{i-D3}9 


xAqi'^ 


63 


no 





97^7(3 £4}? 


x'^Aqc 


4j 


no 




97^7{i-D3L'4}9 






no 





qi^f^zAq 


x6qc 


sr 


no 


^0 


qj^cr3{4Di}q 


di 




no** 





qi^l[-iDi]q 


d2 


8r 


no** 


^0 


qi^l\iD\s\Di}q 



Table 2 

Perturbative renormalization constants 
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13=6.0 


/3=5.6 




8 
3 


-3.16486 


40 
9 


0.989 


0.988 




8 

2^ 


-1.88259 


40 


0.978 


0.977 




-19.57184 


ir 
9 


1.102 


1.110 


x^g 


30 


-35.35192 


2216 
225 


1.215 


1.231 







15.79628 





0.867 


0.857 




8 
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-4.09933 


40 
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0.997 


0.997 
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& 
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-4.09933 


40 

s? 
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0.997 


0.997 


a;^A(7 


-19.56159 


1.102 


1.110 


5g 


1 


16.01808 


-1 


0.856 


0.846 


xSq 


3 


-4.47754 


-5 


0.996 


0.995 


di 





0.36500 





0.997 


0.997 


d2 


7 

6 


-15.67745 


35 

18 


1.116 


1.124 



The continuum operators defined above are ap- 
proximated on a discrete cartesian lattice us- 
ing representations of the hypercubic group that 
eliminate operator mixing where possible and 
minimize the number of non-zero components of 



the nucleon momentum. The operators we have 
used are shown in Table |^, where we have indi- 
cated whether the spatial momentum components 
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Figure 2. Comparison of chiral extrapolations of 
full and quenched calculations of (xq) and {xAq) 
showing agreement within statistical errors 



are non-zero and whether mixing occurs. (Note, 
no* indicates a case in which mixing could exist in 
general but vanishes perturbatively for Wilson or 
overlap fermions and no** indicates perturbative 
mixing with lower dimension operators for Wilson 
fermions but no mixing for overlap fermions.) 

The perturbative renormalization coefhcients 
we have calculated and used in this work are 
tabulated in Table |2| |4| . The factor to convert 
lattice results to the continuum M S sch eme is 
Z{gl = 6/13) = 1 - - B*^^). 



5. RESULTS 

The moments listed in Table 1 were calculated 
|l] on 16^ X 32 lattices for Wilson fermions in 
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Figure 3. Comparison of chiral extrapolations of 
full and cooled calculations of {xAq) showing the 
extent to which instantons reproduce the full re- 
sult 



full QCD at /3 = 5.6 using 200 SESAM config- 
urations at each of 4 k's and at /3 — 5.5 using 
100 SCRI configurations at 3 k's. They were also 
calculated with two sets of 100 full QCD con- 
figurations cooled with 50 cooling steps and in 
quenched QCD at (3 — 6.0 using 200 configura- 
tions at each of 3 k's. Typical chiral extrapola- 
tions for operators calculated with nucleon mo- 
mentum equal to zero are shown in Figure ^ for 
full and quenched calculations of (xq) and (xAq) , 
showing agreement within statistical errors. To 
avoid finite volume errors at the lightest quark 
mass, the SESAM g results were extrapolated 
using the three heaviest quark masses. Table 3 
shows a major result of our work, that there is 
complete agreement within statistics between full 
and quenched results. Statistics with the SCRI 
configurations j|] are not yet adequate to present 
extrapolations in the coupling constant. 

Typical chiral extrapolations for cooled con- 
figurations are compared with the correspond- 
ing uncooled full QCD calculations in Figure ||. 
This qualitative agreement between cooled and 
uncooled results occurs at light quark mass for 
all the twist-2 matrix elements we calculated and 
demonstrates the degree to which the instanton 
content of the configurations and their associated 
zero modes dominate light hadron structure |^]. 
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Table 3 

Comparison of our full QCD and quenched results with other lattice calculations and phenomenology at 
4 GeV in WS 





QCDSF 


QCDSF 
(a = 0) 


Wuppertal 


Quenched 


Full QCD 

(3 pts) 


Phenomenology 

{qval) 


XUc 


0.452(26) 






0.454(29) 


0.459(29) 


0.284 


xdc 


0.189(12) 






0.203(14) 


0.190(17) 


0.104 


XUc dc 


0.263(17) 






0.251(18) 


0.269(23) 


0.180 


X^Uc 


0.104(20) 






0.119(61) 


0.176(63) 


0.083 


X dc 


0.037(10) 






0.029(32) 


0.0314(303) 


0.025 


X'^Uc 


0.022(11) 






0.037(36) 


0.0685(392) 


0.032 


x^dc 


-0.001(7) 






0.009(18) 


-0.00989(1529) 


0.008 


Auc 


0.830(70) 


0.889(29) 


0.816(20) 


0.888(80) 


0.860(69) 


0.918 


Adc 


-0.244(22) 


-0.236(27) 


-0.237(9) 


-0.241(58) 


-0.171(43) 


-0.339 


Auc — Adc 


1.074(90) 


1.14(3) 


1.053(27) 


1.129(98) 


1.031(81) 


1.257 


xAuc 


0.198(8) 






0.215(25) 


0.242(22) 


0.150 


xAdc 


-0.048(3) 






-0.054(16) 


-0.0290(129) 


-0.055 


x^Auc 


0.087(14) 






0.027(60) 


0.116(42) 


0.050 


x'^Adc 


-0.025(6) 






-0.003(25) 


0.00142(2515) 


0.016 


SUc 


0.93(3) 


0.980(30) 




1.01(8) 


0.963(59) 




Sdc 


-0.20(2) 


-0.234(17) 




-0.20(5) 


-0.202(36) 






-0.206(18) 






-0.233(86) 


-0.228(81) 




di 


-0.035(6) 






0.040(31) 


0.0765(310) 





Comparison in Table 3 of our quenched results 
with those of the QCDSF collaboration [|| and 
the full QCD results with those of Wuppertal 
Q shows complete consistency. Since the phe- 
nomenological quantity q^ai |f| does not corre- 
spond precisely to the connected diagrams we cal- 
culate, the most meaningful comparison is with 
u — d differences. The two most physically signif- 
icant discrepancies arising from this table are the 
fact that the difference in first moment is ~ 0.25 
- 0.29 on the lattice and 0.18 experimentally and 
the axial charge is ~ 1.0 - 1.1 on the lattice and 
1.26 experimentally. We have clearly shown that 
these discrepancies do not arise from quenching 
and believe both indicate inadequate treatment 
of the pion cloud of the nucleon due to the small 
physical volume and heavy quark mass. 
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